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Abstract 

We study the dynamics of 4 kinks perturbed by an ac force, both with 
and without damping. We address this issue by using a collective coordinate 
theory, which allows us to reduce the problem to the dynamics of the kink 
center and width. We carry out a careful analysis of the corresponding or- 
dinary differential equations, of Mathieu type in the undamped case, finding 
and characterizing the resonant frequencies and the regions of existence of 
resonant solutions. We verify the accuracy of our predictions by numerical 
simulation of the full partial differential equation, showing that the collective 
coordinate prediction is very accurate. Numerical simulations for the damped 
case establish that the strongest resonance is the one at half the frequency of 
the internal mode of the kink. In the conclusion we discuss on the possible 
relevance of our results for other systems, especially the sine-Gordon equation. 
We also obtain additional results regarding the equivalence between different 
collective coordinate methods applied to this problem. 

PACS numbers: 05.45.Yv, 02.30.Jr, 03.50.-z, 63.20.Pw 
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I. INTRODUCTION 



One century and a half after their discovery, solitons and solitary waves have proven 
themselves ubiquitous in nature, arising in very many physical applications and leading to 
very important advances in applied mathematics [[[],§. Generally speaking, it is often the 
case that the properties of solitary waves are known for certain equations, perhaps integrable, 
that relate to an oversimplified description of different physical systems; subsequently, one is 
interested to learn how these properties are modified if terms initially neglected are included 
as perturbations of the original equation. In most cases, this is a very complicated problem, 
and shedding light on it usually requires the use of approximate analytical approaches. One 
of the most succesful and widely applicable of these approaches is the collective coordinate 
technique or, rather, the family of collective coordinate techniques [§],f|]. The main merit 
of these procedures is the drastic reduction of the number of degrees of freedom involved 
in the problem, from the infinity of them in the original partial differential equation to the 
dynamics of a few degrees of freedom, governed by ordinary differential equations. Quite 
commonly, the reduction is done to a single degree of freedom, which, in general, can be 
identified with the center of the wave X(t) or its velocity V(t), as was first proposed in the 
mid-seventies . This amounts to mapping the motion of solitons or solitary waves to the 
motion of a pointlike (perhaps relativistic) particle with an effective mass [0. Surprisingly, 
this dramatic simplification leads to excellent results for very many systems |3|,[|] . However, 
there are perturbations which, when acting on solitary waves, change not only the position of 
the center of mass of the kink, but also its width, as suggested in , or even some radiation 
can appear in the waves, a phenomenon that can yield the whole collective coordinate idea 



useless as more and more degrees of freedom are excited PJItfl. 

As we have mentioned above, in 1983 Rice || developed a new perturbative method, 
which he applied to two well-known nonlinear Klein-Gordon problems, the 4 and the sine- 
Gordon equations |L[^| , in order to account for variations of the width of their kink solutions 
under perturbations. In Rice's approach, the collective coordinates are the kink center X(t) 
and its width l(t). His results point out that when, in those systems, the kink is subject to 
"some perturbations" (sic) the simple translational motion of the kink center can be coupled 
to an oscillatory motion of the width of the wave, whose frequency he obtained by means 
of a variational approach. Interestingly, for the 4 equation this so-called Rice frequency, 
Qr, practically coincides with the frequency Qi of the kink internal mode (one of the modes 
of linear excitations around the kink, corresponding to a nonzero eigenvalue in the discrete 



spectrum, see, e.g., [PI ). On the contrary, for the sine-Gordon (sG) case, Qr turns out to be 
within the phonon spectrum and, furthermore, the sG kink does not have an internal mode 
(its only eigenvalue in the discrete spectrum is zero, corresponding to a Goldstone mode 
H). For this reason, the deformation of the kink width due to the internal mode has been 
studied extensively in the 4 equation and, among other interesting results, we now know 
that the internal mode is able to store and transfer the energy in 4 kink-antikink collisions 
|T2| |, in the interaction of (ft 4 kinks with impurities JT3 | , or in the case when the 4 kink is 



subject to a periodic spatially modulated potential This exchange of energy between 
the internal and the translational modes or among the internal mode and the modes of the 
impurities explains the resonances that take place in the above perturbed <ft 4 systems. 

In view of these results, a question that naturally arises is the possible existence of 
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resonances of the 4 kink when perturbed by external ac forces, a problem that has not 
been considered so far and that relates to a number of physical contexts. Very recently, 
we have shown |Tj| that a strong resonance arises when the 4 system is subject to an 



external ac force of the form f(t) = esm(5t + 5q) with S close to f2j/2, while in the case 
that 5 = Qi this resonance is weak and even disappears for an appropriate choice of the 
initial kink velocity and the parameters of the driving force. However, the work reported 
|T5f contained mostly numerical results and only an intuitive explanation of this striking 



m 



phenomenon in terms of the collective coordinate equations, and therefore our aim in this 
paper is to provide a full analytical treatment of those equations (for the undamped case) 
in order to better understand the anomalous resonance phenomenon. In addition, in doing 
so we will be already carrying out the same analysis for the sG equation, which may be 
of interest in order to clarify the question of the existence of an internal quasi-mode for 



this equation [ 16| . We deal with these issues along the paper according to the following 



scheme: In the next Section we use the so called generalized travelling wave Ansatz (GTWA) 



17j to obtain the equations governing the dynamics of the kink center and width (which 
is associated to the excitation of the internal mode). In Sec. Ill we thoroughly analyze 
those equations in the absence of damping, identifying all the possible resonances and their 
locations as a function of the equation parameters. We then focus on the most interesting 
resonances, namely those at fli/2 and For the damped case we have not been able to 
solve analytically the equation for l(t), but in Sec. IV we present numerical simulations of 
the full partial differential equation confirming that the resonance at Qi/2 remains for weak 
damping, while the resonance at 5 ~ Qi disappears for any damping value. Additionally, 
in that section we have compared all our analytical results to the numerical simulations, 
finding an excellent agreement. Finally, in the last section we summarize our results and 
discuss their implications for other systems, focusing especially on the sG kink dynamics. 
In an appendix, we prove the equivalence of the GTWA to a different procedure to obtain 
collective coordinate equations, which uses the momentum and the energy of the system 
similarly to the classical approach in ||. 

II. COLLECTIVE COORDINATE APPROACH 

The dynamics of 4 kinks subject to a periodic force f(t) = esm(St + 5o) is governed by 
the equation 

<ptt-<p xx -<P + <f = -f3<Pt + f{t), (1) 

where (3 is the damping coefficient, and e, 5 and <5 represent the amplitude, the frequency 
and the phase of the periodic force, respectively. 

In order to apply the GTWA, first proposed in |T7|], we rewrite Eq. ([]]) in the following 

way: 



6H_ 
SH 



(2) 

if> = -^r + F(x,t,<j>,<lH,...), (3) 



3 



where ip = <ft an d the dot represents derivative with respect to time, F(x,t,<f),<j)t, ■■■) = 
—/3<}> + f(t), and H is the Hamiltonian of the system with the corresponding potential, 
U(4>) = ((p 2 — l) 2 /4, when e and (3 are zero, 

/+oo ,-11 
_ dx{-^ + -<Pl + U(<P)}. (4) 

We now assume that the solution of has the form 

<f>(x,t) = <j>[x-X{t),l(t)], (5) 

and, hence, from the definition of ip we have that 

ip(x,t) =tj) [x-X(t),l(t),X,i\. (6) 

Here <fi describes the soliton shape, whose center will be given by X(t), and where we 
introduced a second collective variable, l(t), which will represent the kink width as we will 
see below. This Ansatz takes into account that I can be different from the Lorentz-contracted 
width due to the action of the external force. Since the internal mode is related with the 
kink width, we can expect that by using this approach we will be able to explain the kink 
motion when the internal mode is excited. However, the linearized problem around the 
initial kink solution tells us that perturbations of the kink can not only shift its position or 



change its width, but also make it emit radiation ||11||. We have not considered this effect 



and therefore the Ansatz is only valid when the external force is so small that practically no 



radiation is emitted. Since the frequencies of the discrete internal mode f2j = J 3/2 and the 



continuum phonon band u p = y 2 + k 2 /2 are apart from each other one can expect that the 
two kinds of modes are excited at different values of the parameters of the external force, 
and then when resonances appear as a consequence of the excitation of the internal mode 
one can in principle neglect the effect of radiation. In any event all our analytical results 
will have to be confirmed later by numerical simulations, which will show whether or not 
this assumption is correct. 



To obtain the equations for our collective coordinates X(t) and l(t) we will follow [17 
(see the appendix for another possible derivation). First, we insert (|§)-(|§) into (0) and (^|), 
and then we multiply the first obtained equation by dip/dX and the second one by d<f)/dX; 
subtracting both expressions and integrating we arrive at the following system of equations: 



r+oo 



r+oo 



I — oo 



dxF(x, t, 



(7) 



Hi 



'" } dx 



with 



^stat 



OC 



d(j) dtp dcpdip , s 

[ ^ ] = dx sT " m 9X' (8) 



5(j> dX 5ip dX J J-oo dX dX 
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where E represents the energy of the system, 7i is the Hamiltonian density and E is 
the static force due an the external field or other solitons, equal to zero for the above 
Hamiltonian. In order to obtain the second equation of motion we proceed as in the previous 
equation and begin by inserting @ and (|) into (0) and @, then multiplying the first and 
second equation by dip/dl and dcf)/dl respectively, and finally taking the difference and 
integrating, thus finding 



+00 



dx [ip,(f)] X + 



di ax 



dl dl 



dx F(x,t 



5 V> Vt, 



(10) 



dl 



where 



K int (l,l,X) 



, &H 

dx ~af 

00 dl 



dE 
'~ffl' 



(11) 



By this procedure, we have obtained two coupled second-order ordinary differential equations 
for X(t) and l(t), where up to now we have not imposed any condition on the soliton 
shape; however, to solve Eqs. flTD-(|iTD we need the explicit functional dependence of (p(x,t). 
Following Rice M we assume that 



(x,t) = <f) [x -X(t),l{t)] =tanh 



x - X(t) 



(12) 



where <po = = tanh[x — Xq/Iq] is the static kink solution of the 4 system, centered at Xo 
and of width l = \J2. Substituting (|12D into (f^)-([TT|) and integrating over x we obtain 



where 



X XI X 

M / - - M l — = F stat (X) - (3M l 0J + F e . 



aMo/07 + Mvk — 



r- 

X 2 
P 



K mt (lJ,X) -fiaMoloj + K 



(13) 
(14) 



dx f(t) 



OX 



-?/(*), 



Tistat 



(15) 



K 



+00 



dx f(t) 



dl 



0, K 



int 



dE 
~dl' 



(16) 



*4t^W + HK)' (l7) 

with a = (tt 2 - 6)/12, q = 2 and M = 4/(3/ ). Denoting P(t) = M l X /I, the Eq. (p|) can 
be written as 
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It is interesting to note that this equation may be obtained as well by applying the McLaugh- 
lin and Scott procedure || with one collective variable only corresponding to the center of 
the kink. As was shown in |Tj|, for the sG kink dynamics this is an excellent description 
of the kink motion, and we have verified that it also describes the 4 kink motion under ac 
forces away from the resonances we will find and discuss below. Its solution is given by 



q e 



P{t) ((3 2 + 5 2 ) 
+ exp(-pt)<P(0) 



[5 cos(5t + 5 ) - (3 sin(5t + 5 Q )] + 
q e 



(P 2 + 5 2 



[(3sm(5 ) - 5 cos(5 )] >■ 



(19) 



From Eq. (0) the equation that holds for the width l(t) is 



a 



I 2 - 211 - 2(311 



I 2 
'o 



1 + 



M 2 



- 1. 



(20) 



Note that the term P(t) 2 in Eq. fl2D| ) involves two frequencies 5 and 28 when f3 = 0; whereas 
for (3 and after some transient time the only frequency that remains is 25. Furthermore, 
the term with frequency 5 vanishes for (3 = and an appropriate choice of initial parameters. 
This equation represents a nonlinear, damped and parametrically excited oscillator. To solve 
Eq. (|20|) we provide the following change of variables, proposed in [f[£(| , l(t) = g 2 (t), which 
transforms the above equation into an Ermakov-type equation (or Pinney-type) [see [OIJ and 
references therein]. Then, the equation for g(t) reads 



9 + 13 9 
9(0) = 



n\ 2 

2) + 



n 



2M 



P 2 



1 



h + 0, g(0) 



2Vfs' 



4ag 3 ' 



(21) 



1 / y/al s in the case when 
1.2247. We have not 



where Q = 1/y/alo = 1.2452 is equal to the Rice frequency Qr 
the kink initially is at rest; this agrees within 1.7% with Oj = y3/2 
been able to solve analytically Eq. (0), except when (3 — 0. The next section is devoted to 
a detailed analysis of that case; we will come back to the nonzero (3 problem when discussing 
our numerical results in Sec. IV. 



III. UNDAMPED KINK: = 



When = 0, Eq. (gT|) reads 



2. 



n 

2M Q 



P' 



Aag 3 



(22) 



where P(t) is given by P(t) = A + qecos(5t + 5 )/5, with A = M 7om(0)// s — qecos(5 )/S. 
We thus see that the function P(t) 2 in Eq. (|22D involves trigonometric functions with fre- 
quencies 5 and 25 if and only if A ^ 0. On the contrary, when A = the only frequency that 
remains in the function P(t) 2 is 25. Interestingly, we note that the relation A = coincides 
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with the condition for the oscillatory motion of the center of the kink, obtained by using the 
McLaughlin and Scott approach in the absence of dissipation ]n| . The solution of Eq. ( |22|) 
El is 



»<*> = (23) 

where V\(t) and v 2 (t) are two independent solutions of the linear part of Eq. (|22|), W = 
v\v 2 — v 2 Vx is the Wronskian, which in this case is a constant and can be calculated by the 
initial conditions for Vi (i = 1,2) which are t>i(0) = \fT s , V\(0) = 1(0) / (2^ l /%) , v 2 (0) = and 
^2(0) = const ^ 0. 



A. A = 0: Resonance at 5 « Or/2 

If one denotes t = 5t + 5q and sets A = in the linear part of Eq. fl22"|), after some 
manipulations we arrive at the following Mathieu equation for the V{ functions: 

v"+ [a + 20 cos(2r)]v< = 0, 

n\ 2 q 2 e 2 (24) 







g 2 e 2 " 




1 H 


25 2 M 2 _ 



2 ; 



257 45 2 M 2 

where prime denotes the derivative with respect to r. Notice that the initial conditions for 
Vi(r) become Vx(8 ) = y/T s , v£(<f ) = Z(0)/(25 v / Q, v 2 (5 ) = and v' 2 (5 ) = v(0)/5. The 
solution of Eq. ( |2"4"D (see P^| ) for V\(r) and v 2 (r) can be expressed as a linear superposition 
of the two Mathieu functions ce u and se u with a non- integer index v, i.e., 



where 



and 



Ui(r) = A iC e u (r,-9) + B iS e u (r, -9) , 1 = 1,2, (25) 



4 = %, -Si = (26) 



A = ce„ (S , -9) se' u (5 , -9) - ce' u (8 , -9) se u (5 , -9), 
A Al = Vi(S ) se' u (S , -9) - v'^Sq) se u (5 , -9), 
A Bi = v'^80) ce u (S , -9) - Vi(5 ) ce' u (S , -9), 
with the constraint (characteristic curve for Mathieu functions) 

a = , 2 + ^-i- Ty 2 + O(^). (27) 

From Eqs. (|23|), ( |25|) and (^6|), and taking into account that r = 5t + <5 we obtain that the 
kink width l(t) is given by 

l(t)=9 2 = vl(t) + I ^ r y 2 (t), (28) 



7 



Vi(t) = At ce u (St + 5 , -9) + Bi se v (St + S ,-9), i = 1, 2; 



(29) 



where 

Vi(0) se u (6 Q , -9) - Vi(0)se u (S ,-9) 



A 



[ce u (S , -9) se u (S 0} -9) - ce v (6 Q , -9) se u (5 , -9)] ' 
B _ Vi(0)ce v (5 ,-6) - ^(0)^(50,-9) ^ 

[ce„ (S , -9) se u (S , -9) - ce v (S , -9) se u (S , -9)] ' 



W = -V^ 2 (0), (31) 
and the characteristic curve, Eq. ( P7|) , for our initial parameters can be written up to order 



e 2 as 



^_gWgo) 
2v 2M 2 1 2 n R 



To obtain a better approximation, we need to take into account more terms of the above 
series, so in Fig. [I] we have plotted, as solid lines, the characteristic curves obtained numer- 
ically with Mathematica 3.0 [p3|l , for v = 1/2 and v = 3/2 when u(0) = and So = tt/2. 
Notice that when v = m + p/s is rational, with m an integer number and p/s a rational 
fraction (0 < p/s < 1), Vi(t) or v 2 (t) are 2ns periodic functions, if p is odd, and 7rs-periodic 
functions, if p is even; whereas for irrational v both functions will be non-periodic, but 
bounded solutions p2| . 

If we know l(t), then from the solution P(t) [Eq. dl9|)1, we can calculate the velocity 
X = P(t)l(t) I '(MqIo) for the kink center. Since the momentum P(t) is a periodic function 
and the kink width l(t) at least is a bounded function, X is bounded as well. For instance, 
if we take v = 1/2 in fl3"2]), the frequency of the external force 



g 2 cos(26 ) 2 4 

5= ^-i^r +0(e) ' 



is very close to fl^ ~ fli if e 1 (see Fig. [I]), and /(t) is a function of the square of the 
47r/<5-periodic solutions sei/ 2 and cei/ 2 , respectively, so the velocity of the kink center, X(t), 
and the energy E(t) will also be bounded functions for S = S(e). 

If we try to find, e.g., 27r-periodic solutions of (|24"D we obtain that in this case the two 
independent solutions are related to the integer Mathieu functions sei and cei. However, 
these two solutions appear when the characteristic curve a(9) is 

„ 9 2 9 3 9 A 

a=l + 9 h 

8 64 1536 

for the even Mathieu function cei and 

„ 9 2 9 3 9 A 
a = 1 — 9 



64 1536 



for the odd Mathieu function sei. Since the values of a are different for each of the function 
sei and cei, these functions are not solutions of the same equation except when 9 = 0. 



S 



Although sei and cei are not solutions of Eq. (p4j ) [recall that in order to find l(t) we need to 



calculate two independent solutions of Eq. (24)], the characteristic curves of these solutions 
separate the unstable and stable solutions of the Mathieu equation (p4|), and can be rewritten 



as 



and 



Or 
2 



5 

+ 2 



q 2 cos(25 



+ 



(33) 



q 2 cos(25 ) 
2M 2 1 2 Q R ~~ 4M 2 -f 2 Q R 



(34) 



respectively. Interestingly, we note that the width of the unbounded region, 

A5 = 5 + - <L 



q 2 e 2 

2„,2f 



2M 2 7o 2 Q B ' 

decreases when the initial velocity increases, and consequently we will focus on the case 
u(0) = in our numerical simulations below, as the resonance is then easier to observe. 
The curves 5 + and <5_ [Eqs. (|33"D and ([341) 1 are plotted in Fig. [I] for zero initial velocity 
and 5q = 7r/2, where the shadowed region represents the region where l(t) is unbounded. 
Analogously, we can obtain other characteristic curves, related either with integer Mathieu 
functions se n and ce n (n e N), or with non-integer Mathieu functions se u and ce u . In 
view of this, the above analytical results lead us to expect that if 5 is close to Q R ~ Qf, 
l(t), u(t) and the energy should be oscillatory (or, at least, bounded) functions, whereas 
if 5 ~ Qr/2 « Qj/2 the kink width l(t) increases indefinitely and since the energy and the 
velocity are proportional to l(t) [see Eqs. (|T9|) and (|T7|) 1 we should observe that these 
functions increase with time as well. 



B. A / : Resonances at S « Qr/2, £l R 



When A 7^ 0, the linear part of Eq. ( |2"2"D becomes, in a manner similar to the procedure 
leading to Eq. (E4l), a more general Mathieu equation ||24||, namely 



b + 2B x e 2 cos(2r) + 29 2 e 2 cos(4r) g = 



(35) 



where 



Q\ 2 
6J 

Xq 2 



q 2 e 2 

+ 25 2 M 2 + i r2 



Ml 



5 2 Ml 



0* 



A5 2 M^ 



and the prime denotes the derivative with respect to r = (8t + S Q )/2. According to Floquet 
theory | 25| , Eq. (|35|) has normal solutions of the form g = exp(crr) ( l ) (r), where $(t) is 
a 7r-periodic function and a is the characteristic exponent. Expanding $(t) in a Fourier 
series, the function g can be rewritten as g = Y,t=-oo ( ^ > n ex P[(°" + 2m)r], where <3> n are the 
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coefficients of the above series. Therefore, the transition curves separating stability from 
instability correspond to a = (7r-periodic solutions) and a = i (27r-periodic solutions) [26 



Now we will apply the method of the strained parameters [53] to determine the values of 
b, d\ and 9 2 corresponding to these values of the characteristic exponent. First we assume 
the solutions of fl55|), which have period tt or 2n and the transition curves b = 6(e), can be 
written in the form of perturbation expansions as 



Second, substituting 
of e we obtain 



g(r)=g + e 2 g 1 + e A g 2 + O(e (i ) 
b = b + e 2 b 1 + e% + O{e 6 ). 



(36) 



into (|35j) and equating to zero the coefficients of the same order 



£go = 0, 



(37) 



Cg x = -hgo - 26»i cos(2r)5f - 2# 2 cos(4r)sf , 



(38) 



£#2 = -hgi - b 2 g - 29 \ cos(2r)5fi - 29 2 cos(4r)gi 



(39) 



where £ represent the second-order linear differential operator C = — j + 6 . The solution 



of Eq. (§7|) is 



g = A cos y b r + B sin I y b r 



(40) 



where bo = 4n 2 for the 7r-periodic solutions and bo = (2n — l) 2 for the 27r-periodic solutions, 
n is an integer number and A and B are constants. 

For the 7r-periodic solutions we will analyze the case n = 1, corresponding to bo = 4; we 
note that the case n = is not possible due to the definition of our parameter b. When 
bo — 4, g — v4cos(2r) + £>sin(2r). For g\ to be a 7r-periodic function is necessary to 
eliminate the secular terms in Eq. (|3lf) , imposing either b\ = 9 2 and A = or b\ = —9 2 and 
B = 0. Then, the transition curves corresponding to the solutions 



g + = Bsm(2T) + Be 2 



^ sin(4r) + ^| sin(6r) 



(41) 



and 



g _ = A cos(2r) + At 2 
can be shown to be 
or, alternatively, 



"^i i ^ (a \ i ^ 2 ^« ^ 
cos(4r) H cos(6r 

4 12 V ; 32 V ; 



6 = 4 ± e 2 ^ 9 



(42) 
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5, 



Or 
2 

2 



qu(0) cos(5 ) 

Mo7o 
qu(0) cos(5o) 

M)7o 



3g 2 



e + 



4M 2 7o 2 n R 



+ 



+ 



g 2 cos(25 ) 
q 2 cos(25 ) 

2^,2( 



(43) 



4M 2 7o 2 ^ 2M 2 7o 2 O ii 



respectively. Since these curves start from e = and b = 4, we are analyzing the transition 
curves when the frequency of the ac force is close to half of Qr. Hence, the resonance at 
5 ~ Qi/2 ~ appears much as in the case A = even when A ^ 0. Let us remark that, 

for given e, the distance between these two curves 5 = 5(e) is 



g 2 e 2 



2M 2 7o ^V 

and when the velocity increases, the unstable region is more narrow, as we found above in 
the case when A = 0. 

The other transition curves are related to the 27r-periodic solution of Eq. (|3~TD, go = 
Acos(t) + -Bsin(r), which appear when bo = 1 (n = 1). In this case g\ will be a periodic 
function if b\ = Q\ and A = or b\ = —9i and B = 0. Hence, the transition curves starting 
from 6=1 are 

b= l±e 2 #i, 



or, equivalently, 

qu(0) cosQq) 
d T = ll R — e + 

M 7o 

and correspond to the solutions 



q 2 cos(25 ) 
2M 2 1 2 n R ' 4M 2 7o 2 ^ 



+ 



e 2 T 



q 2 \ 



2M 2 ^n R 



(44) 



BsmM + Be 2 



9i 



@2 @2 

— sin(3r) + — sin(5rj 



(45) 



y4cos(r) + Ae 



— - — cos(3r) + — cos(5r) 



(46) 



respectively. These solutions and its corresponding characteristic curves are related with 
the driving frequency 5 close to fl R . Finally, in this case we have 



A5 = 5 + - 6- 



q 2 Xe 2 



which means as above that for fixed A, AS decreases when u(0) is increased. 

In summary, we have found here that when A ^ and we drive the system with a 
frequency close to fl R the solutions will be unstable in the region between the curves 5 + and 
5-. Notice that for A = we recover the previous results, i.e., there are not any resonances 
at 5 ~ Q R . 
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IV. NUMERICAL VERIFICATION 



The results we have obtained in the previous section have been derived within the collec- 
tive coordinate assumption that we can describe all the kink dynamics by the two variables 
X(t) and l(t), all other degrees of freedom being negligible. As there is no way to know a 
priori that this is indeed the case, we have to verify this hypothesis by numerical simula- 
tions of the full partial differential equation. Furthermore, we have not been able to solve 
the collective coordinate equations for the damped (/3 7^ 0) case; it is reasonable to expect 
that in this situation phenomena similar to those found for the undamped case will arise, 
but in so far this is not checked this assertion remains a conjecture. 

In view of the above considerations, we have computed the numerical solution of the 
partial differential equations (0) and ([5J) by using the Strauss- Vazquez scheme p^] and 
choosing a total length for our numerical system of L = 400, with steps At = 0.01, Ax = 0.1, 
free boundary conditions, and simulating up to a final time equal to 25 000 with a kink at 
rest as initial condition. In addition, we have fixed the amplitude e = 0.01 and the phase 
<5o = tt/2 of the ac force and then we have changed the value of the driving frequency 5 
in order to see the resonances. We note that there are many other parameters we could 
change, such as the initial velocity, the driving amplitude, or the driving phase, but as our 
main goal in this section is to assert the validity of the general analytical results obtained 
above we prefer to concentrate on a few mentioned. In all our simulations we have 

chosen the initial parameters in such a way that A = 0, because when A 7^ the kink moves 
to the right or to the left and then for large times the kink will leave our finite system. For 
the aforesaid values of initial parameters the width of the resonance region AS = 5 + — 5- 
predicted at 5 « Vt R /2 ?s fJj/2 in the section III (A) [see Eqs. (p3|) and (|34])1, is of the order 
of 10 -4 . For this reason we have explored the regions around fli and fli/2 in an interval 
of that order. Finally, we want to mention that with the Strauss- Vazquez method one can 
compute very accurately the position and the velocity of the kink center using the integrals 
of energy and momentum [28|] , so it is a good numerical method in order to compare to our 
analytical predictions. 

First of all, in order to verify the results obtained by means of the GTWA with two 
independent collective coordinates, the position and the width of the kink, we have studied 
the region around 5 ~ Qi = 1.2247, i.e., the driving frequency for which we do not expect 
any resonances when A = 0; hence, the width and the velocity should be, at least, bounded 
functions. We have verified that this prediction is in excellent agreement with the numerical 
simulations, an example of which is shown in Fig. |2|. In this plot, we can see that the velocity 
is indeed an oscillatory function and even for large times our theoretical approach describes 
correctly the evolution of the velocity of the center of mass of the kink. Other values close 
to 5 = 1.21 behave very much like the presented example. 

The next test of our analytical results has to be, of course, the existence of a resonance 
near Q.i/2. For frequencies 5 around Qr/2, our collective variable approach has predicted 
that the width l(t) increases unboundedly and hence the velocity and the energy should 
increase as well. Again, the prediction is fulfilled, but for fli/2 instead of VLr/2: In the 
numerical simulations we have observed that in this case a resonance takes place when 
5 ~ fii/2 = 0.6124, As a specific example, the results for the kink energy for 5 = 0.6104 
and 5 = 0.608 are plotted in Fig. |3], clearly showing a resonant increase of the energy in 
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the former case and oscillatory behavior in the latter. Interestingly, at the resonance the 
velocity also does not have an oscillatory behaviour as depicted in Fig. |], in contrast to 
the behavior off-resonance we found in Fig. |2| for 5 = 1.21. Such a non-periodic evolution 
of the kink velocity implies that the kink motion is chaotic as we first found in |15J and, 
eventually, the kink begins emitting radiation. As we discussed in [IS], we believe that this 
phenomenon has its origin in the energy transfer from the kink internal mode to the rest of 
modes in the system, a mechanism already demonstrated in [I21.H]; intuitively, the reason 



for that is that the collective coordinate prediction that the kink width grows without limit 
can not be physically true due to the own nature of the kink, and hence, when the internal 
mode excitation reaches large values the energy ends up being transferred to the rest of the 
available modes (translation and radiation). 

After checking our analytical results for the undamped case, we now have to turn to the 
damped {(3 ^ 0) dynamics, in order to find out whether the same phenomena arise there. 
Simulations for the damped case show that the energy of the system also increases when 
the resonance takes place (see Fig. |5j), but in this case the energy is bounded due to the 
dissipation, reaches an asymptotic value, and thus does not increase indefinitely as in the 
absence of damping. For comparison of both situations we have computed the mean energy 
in the time interval 10 000 < t < 25 000 with and without damping, and our results are 
summarized in Figs. |6| and 0. For (3 = 0.001 (see points on the solid line in Fig. []) we have 
found that the energy increases at 5 = 0.6103, whereas for (3 = (see the points on the 
dashed line in Fig. the resonance frequency is 5 ~ 0.6102. 

As mentioned above, the numerical solutions of the full partial differential equation show 
that the resonance at S ~ flj/2 also takes place for (3^0, whereas when S ~ Q{ (see Fig. |7|) 
the resonance disappears. Let us recall that in all our simulations we have chosen u(0) = 
(and hence the resonant region should be the widest possible one) and 5 = it/2, in such 
a way that A = 0. In the case when A vanishes we can not expect any resonances at 
5 ~ Qi. Nevertheless, in Fig. [7] it is clear that the energy increases weakly in this region. 
We believe that this maximum of the energy is probably due to the fact that the condition 
for the suppression of the resonance at 5 ~ f2j is extremely difficult to fulfill numerically. 
On the other hand, it is also possible that this condition, which has been obtained within 
the collective coordinate framework, is only approximately true, and therefore when we 
choose A = we are beginning with an initial condition close to the one needed for complete 
suppression of the resonance but not exactly there, and hence we would still see the small 
bump in Fig. |7|. Another possible reason for this peak is that it is possible that when A = 
at 5 ~ Qi ~ a resonance appears in the next order corrections (see the work of Segur 
||29|| , in which the frequency 2f2j arises, and could then be parametrically excited by 5 = Qi). 
Were this the case, this "weaker resonance" , that completely disappears for (3 ^ [see the 
solid line in Fig. IjJ, is the resonance that we can expect either when A ~ in the partial 
differential equation or for A = and large enough times. 



V. CONCLUSIONS 

We have shown analytically and numerically that the internal mode of the 4 model 
can be excited if we drive the system with an ac force of a frequency close to 1/2 of the 
internal mode frequency. This is a very surprising result as the driving we are applying 
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to the system is not parametric. At resonance, as a consequence of the increment of the 
energy of the system, the kink initially at rest begins to move chaotically [15] and also begins 
to radiate, i.e., the energy is transferred from the internal mode to the translational and 
radiational ones. The chaotic motion is confirmed by noting that, when the internal mode 
is excited at 5 ~ in the discrete Fourier transform (DFT) of the kink's velocity [see 

Fig. U we see some frequencies in the low frequency part of the spectrum aside from the 
frequency of the ac force, 5. What is more important, we have presented a full analysis 
of the collective coordinate theory for this problem, whose validity has been undoubtedly 
confirmed by numerical simulations of the full partial differential equation. In particular, 
in the absence of damping we have been able to find analytically all the resonances as well 
as their dependence on the kink initial conditions. On the other hand, we have also shown 
that the resonance at 5 ~ Qi can be suppressed if we add a small damping to the system 
or by an appropriate choice of initial parameters of the ac force. We have verified in our 
simulations that for 5 ~ fli the lower radiational modes are excited and for this reason the 
energy in Fig. |7| increases when 5 increases. Of course, this excitation of the lowest phonons 
cannot be explained within the present collective coordinate theory; a much more involved 
approach including phonon effects (as in ]30|] for the sG model) would probably account for 
that resonance. 

Beyond the application of our results to the 4 kink dynamics, we believe that the same 
phenomenology will arise in other systems for which internal modes are present as, for exam- 



ple, in the double sG equation ||31|| . The analysis we have presented can be straightforwardly 
extended to this system, and it is quite likely that similar resonance effects will arise, as well 
as in other models with the same feature of internal modes. A related interesting question 
concerns the applicability of these results to the sG equation. As we mentioned in the in- 
troduction, our calculation applies directly to the sG equation, but in this system the kinks 
do not posses internal modes. Therefore, the identification f2j « Qr we have done for the 4> 
kink is not available any more, and, what is worse, the putative Qr lies within the phonon 
band. However, it has been reported that close to Qr sG kinks might exhibit a quasi-mode 
|T6| with a very long life. If such a mode indeed exists, which is still in doubt since no other 
reports of its existence have been published, it could be found by looking at resonances, as 
we have done in this paper, at half its frequency, which would be outside the phonon band. 
Work along these lines is in progress p2| . 
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APPENDIX A: EQUIVALENCE OF COLLECTIVE COORDINATE METHODS 

In this appendix we demostrate that the GTWA, the method we used in section II 
to obtain the collective coordinate equations, is equivalent to using the variation of the 
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momentum and the energy of the system (|T]). First of all, from Eqs. (0) and ( |10D we 
can obtain the variation of the momentum and the energy of the perturbed 4 equation. 
Combining the first three terms of Eq. (|7|) and doing some straightforward transformations 
the equation for the momentum can be recast in the form 



-dt = ~U dx dx + L d ^tAAu-)^ 



(Al) 



+00 



P(t) = - / dxfcfa. 
J — 00 

Furthermore, if instead of Eq. (|T0|) we add Eq. (0) rewritten as 

- t, 0, <f> u ...)|| = 0, 

multiply this equation by X, and the expression obtained from flIU|) by I, we conclude that 
Eq. (p7D|) is equivalent to 



— = / dxF(x,t,<p,qj t , 
at J-00 at 



H(t) = j ^ cfejf + f + (A2) 

This implies that when we apply the GTWA in those systems we are, in fact, varying the 
momentum 

dP d f+°° r+00 

-jr = --77 / dx<f) x <f)t = - dx [4> xt (pt - 4>x4>tt} , (A3) 
dt ctt J— 00 J— 00 

and the energy 

£ = i/>{5* + 5* + ,/ <0}= < A4 > 

dx {<fi x (f) xx + (j> t (j) t t + U'(<j))<j> t } , 



where 

0(a;,t) =#r-X(t),Z(t)], & = iP(x, t) = V[x - X(t), l(t),X(t), l(t)]. (A5) 



Moreover, we can also start from Eqs. ( |A3| ), ( |A4j ) and ( |A5|) and obtain (|7|) and ( P~0| ) 



Doing that, we only need to take into account that there are at least two different ways of 
transforming the above integrals: in one of them one substitutes (|A5[) in (|A3|) and (|A4|) , so 
that 
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dP 
HI 

dH 
~dt 



+00 



dx 



x 



dx 



+ 1 



OXdX 

\dx + dxdx 



ox dl 



(A6) 



dX dl 



I L J dl dX dl dl 



(A7) 



whereas, in the other one, we first substitute the systems of Eqs. in (|A3| ) and (|A4|) 

and then assume that <fi(x,t) and ip(x,t) satisfy (|A5|) . Thus, 



dP 
~dt 



+^ , on 

-00 dx ox + 



dB_ 
~~dt 



X 



d(j) 



dX 



+00 Qfj) 

dxF(x,t,<p,(f) t , ...)—, 

-00 CA 

d(j) 



dx F(x,t,0,0t, ...)-— + / / dx F(x,t,(f),(j) t , ...) — 



dl 



(A8) 



(A9) 



Equating (|A|) and flA"8p , and flAl) and (|A9|) , respectively, we find that (|A3|) and (|A4| ) 
become (0) and fllCf), respectively. 

Let us note that the connection between the GTWA and the the variation of the mo- 
mentum and energy gives a physical interpretation of such a technique, and furthermore, 
this second method leads more directly to formulas for X{t) and l(t). Furthermore, if the 
solution of Eqs. (@)-(||) for F(x, t, </>, fa, ...) = —f3fa-\-f(x,t) is a kink centered at X(t) whose 
width, l(t), depends on time, then Eqs. (|A3|) and (|A4]) become 



and 



P(t)l(t) 



+1 



M l 
~P(t) 



dP 
~dt 



/+00 
dxf(x,t)<f> 
-00 

21 i 2 



2 1 

+ -aM l 



2M l 2 

+00 



+(3a l jl + 



d6f(W + X(t),t) 9 fa 



+ 

1 ^ 



(A10) 



(All) 



/ 2 



+ 



0. 



respectively. The first bracket on the right hand of ( |A11| ) vanishes because of ( |A10| ), so, the 
solution of (|A11|) is / = or 



a[l 2 - 211 - 2(311} 



+ 



r 2 
'0 

2l(tf 



1 + 



Ml 



Mai 



- 1 + 

d9f[ei + x(t),t] e<f> . 



(A12) 



In such a way, ( |A10| ) and ( |A12| ) can be interpreted as the equations which the collective 
coordinates, X(t) [P(t)] and l(t), satisfy, and they are obtained from dP/dt and dH/dt, 
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respectively. For example, if f(x,t) = esm(5t + So) in (|A10|) and (|A12|) , we recover the Eqs. 
(|18|) and (|20D , obtained by applying the GTWA in section II. 

As a final remark, we point out that, since we have not used in any moment the explicit 
form for the potential function U(<p), the equivalence between these two aforesaid methods 
remains true for any other perturbed nonlinear Klein-Gordon equation of the form 
when the Hamiltonian of the system is of the form (S) . 
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FIG. 1. Characteristic curves 5 = 5(e) corresponding to Eqs. (|32| ) [upper solid line], (|33| ) and 
( |34| ) [upper dashed lines] for zero initial kink velocity and 5q = ir/2. The lower solid line is related 
with the solution (f2gj)-(|3l]) when v = 3/2. The lower dashed lines are the characteristic curves of 
the integer Mathieu functions with n = 2. Shadowed regions represent unstable solutions of Eq. 
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FIG. 2. Verification of the collective coordinate method in the absence of damping. The solid 
line has been obtained from Eqs. ( |l8|) and (|20|) and represents the values of the velocity of the kink 
center as given by the collective coordinate approach; the points are the velocities of the center 
of the kink from the numerical simulations of Eqs. @-(|3|) by using the Strauss- Vazquez scheme, 
starting from a kink at rest and with force f(t) = 0.01 sin(1.21i + vr/2). The velocity function 
is represented only from t = 24900 to the final time t = 25 000, but the same behaviour of this 
function is observed during the whole run. 
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FIG. 3. Kink energy when the frequency of the driving force 5 is close to Oj/2. The parameters 
of the numerical simulations of the full partial differential equations, Eqs. are (3 = 0, 

e = 0.01, 5o = tt/2, n(0) = 0, 5 = 0.6104 (upper curve, resonantly increasing energy) and 5 = 0.608 
(lower curve, non-resonant oscillations), L = 400, Ax = 0.1, At = 0.01. Notice that half of the 
Rice frequency is £Ir/2 = 0.6226 and half of the frequency of the internal mode is Q.i/2 = 0.6124. 



0.10 




21 



0.030 



0.025 



O 020 
f 0.015 



E 
< 



0.010 



0.005 



0.000 



0.00 



0.05 



0.10 0.15 
(0/(2 71) 



0.20 



0.25 



FIG. 4. Upper panel: Velocity of the kink center computed from the numerical solution of the 
partial differential equation when the frequency of the driving is close to half of the internal mode 
frequency (<5 = 0.6104), (3 = 0, e = 0.01, 5q = it/2 and u(0) = 0. After some time the velocity 
function departs from its oscillatory behavior, which was transiently exhibited at the beginning of 
the run, and develops chaotic features. Lower panel: Discrete Fourier transform of the signal in 
the upper panel. 
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FIG. 5. Resonance in the presence of damping: Kink energy at frequency 5 = 0.6103 obtained 
from numerical simulations of Eqs. (0) and (0) with e = 0.01, (3 = 0.001, n(0) = and 5q = ir/2. 
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FIG. 6. Resonance at fJj/2. Mean energy of the system, computed averaging the energy in the 
numerical simulations of the partial differential equations for t > 10 000 with final time equal to 
25 000; e = 0.01, Sq = 7r/2, u(0) = 0. Dashed line: undamped case, = 0; solid line: damped case, 
= 0.001. 



1.28 
1.23 



o5 1-18 
0) 

iS 1.13 



CT5 
0) 



1.08 



1.03 



0.98 



1.21 




1.22 



1.23 



1.24 



1.25 



1.26 



FIG. 7. Resonance at Jlj. Mean energy of the system, computed averaging the energy in the 
numerical simulations of the partial differential equations for t > 10 000 with final time equal to 
25 000; e = 0.01, Sq = vr/2, u(0) =0,0 = (dashed line) and = 0.001 (solid line). 
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